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ABSTRACT 

Alms. The statistical equilibrium of neutral and ionised silicon in the solar photosphere is investigated. Line formation is discussed and the 
solar silicon abundance determined. 

Methods. High-resolution solar spectra were used to determine solar log gfssi values by comparison with Si line synthesis based on LTE and 
NLTE level populations. The results will be used in a forthcoming paper for dilferential abundance analyses of metal-poor stars. A detailed 
analysis of silicon line spectra leads to setting up realistic model atoms, which are exposed to interactions in plane-parallel solar atmospheric 
models. The resulting departure coefficients are entered into a line-by-line analysis of the visible and near-infrared solar silicon spectrum. 
Results. The statistical equilibrium of Sii turns out to depend marginally on bound-free interaction processes, both radiative and collisional. 
Bound-bound interaction processes do not play a significant role either, except for hydrogen collisions, which have to be chosen adequately for 
fitting the cores of the near-infrared lines. Except for some near-infrared lines, the NLTE influence on the abundances is weak. 
Conclusions. Taking the deviations from LTE in silicon into account, it is possible to calculate the ionisation equilibrium from neutral and 
ionised lines. The solar abundance based on the experimental /-values of Garz corrected for the Becker et al.'s measurement is 7.52 ± 0.05. 
Combined with an extended line sample with selected NIST /-values, the solar abundance is 7.52 ± 0.06, with a nearly perfect ionisation 
equilibrium of Alog£Q(Sill/Sil) = -0.01. 
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1. Introduction 

Neutral silicon gives rise to only a very few and predominantly 
weak spectral lines in the visible. For that reason it is some- 
times regarded as unspectacular and consequently ignored in 
many abundance analyses. There are at least three important 
facts that put Si back into perspective: 

- Silicon is commonly attributed to the a elements that are 
assumed to be produced in massive type II supernovae. In 
this case the Galactic chemical evolution of Si should dif- 
fer significantly from that of iron resulting in logarithmic 
overabundance ratios [Si/Fe] > 0. 

- Silicon belongs to the most abundant metals. It competes 
with Mg and Fe, which also have similar ionisation ener- 
gies. It is therefore important as an electron donor in rel- 
atively cool turnoff stellar atmospheres, where elements 
such as H, C, N, O and Ne with even higher abundances 
cannot contribute. 



- Due to its condensation behaviour. Si is found in many in- 
terstellar and interplanetary dust particles and meteorites. 
Therefore, Si is the reference element for abundance de- 
terminations from solar system meteorites. This is different 
from the stellar abundance reference, hydrogen. 

All three points cannot hide the problem every astronomer 
encounters when trying to determine the Si abundance in cool 
turnoff stars. Even the line spectrum observed in metal-rich 
stars is poor, which may be part of the reason why the labo- 
ratory work is also far from complete, and some terms and line 
transitions have not yet been identified unambiguously. On the 
other hand, spectra of hotter stars such as those of types B or O 
display a multitude of spectral lines of up to three consecutive 
ionisation stages. Si ii. Si iii, and Si iv. 

Neutral silicon has only a single strong line in the blue part 
of the spectrum, i.e. the 3p 'S^ - 4s 'P° 3905 A line. That line 
appears in a spectral range that is not easy to observe because 
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- cool stars have their maximum flux in the red, 

- spectrograph coatings are optimized for the visual, and 

- CCD detectors have a lower QDE in the blue. 
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Fig. 1. Term energy diagram of the silicon model atom used in this analysis with Sin (top) and Sii (bottom). Sin quartets are 
neglected. Allowed transitions are continuous, the forbidden Si 1 4103 A line is dotted. 



Additionally, the line is blended by a number of other metal 
hues, and line synthesis requires a careful representation of the 
blend components, in particular in metal-rich stellar spectra. 
A second important line is the intercombination line, 3/? 'S,, - 
4i^Pj at 4103 A. Different from other allowed transitions in 
the visible and the near red, this line is relatively strong (com- 
parable to the Mg I line at 4571 A), and at typical S/N of stellar 
spectra the detection becomes difficult at metalhcities [Fe/H] 
< -2.5. 



The situation becomes increasingly worse in metal-poor 
stars, even more so if the stars are at the hot end of the turnoff 
with effective temperatures around 6200 to 6500 K. Here, Si i 
as a typical minority ion becomes nearly fully ionised leaving 
practically no neutral silicon atoms in any of the excited lev- 
els. Since the triplet lines arising from the 3p^P ground state 
are all in the satellite UV, they cannot be observed from the 
ground. The few near-infrared lines are absorbed from levels 
above 5 eV and also become invisible in the spectra of metal- 
poor turnoff stars. 
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Fig. 2. Photoionisation cross-sections of the metastable terms of the Sii model atom as given in the TOPbase. All three cross- 
sections are significantly larger than the corresponding hydrogenic approximations (dotted curves). The absorption edge is located 
at the right boundary of the individual frames. 



Whereas the solar spectrum still provides enough absorp- 
tion lines for a thorough abundances analysis, this is not the 
case for hot metal-poor stars. Even the 4103 A intercombina- 
tion line cannot be used there for an abundance analysis, and 
silicon abundances of such stars will have to rely on a single 
line, /13905 A. This makes an LTE abundance analysis partic- 
ularly prone to systematic errors due to atomic interaction pro- 
cesses that tend towards a net over- or underpopulation of the 
ground state. As we shall see below, the Bp ^P^ j ^ levels and the 
metastable 3;? 'Dj and 3/? 'Sq levels are strongly depopulated 
by photoionisation, and a NLTE analysis of the line formation 
problem is therefore unavoidable. 

As mentioned above, silicon spectra of cool stars were not 
investigated often. Only one NLTE analysis has appeared to 
date. Wedemeyer ('2001') reported the influence of deviations 
from LTE upon the Si abundances in the Sun and in Vega with 
no reference to line formation in metal-poor stars. A number of 
LTE analyses, however, cover the whole range of metal abun- 
dances, from solar to the most extreme metal-poor halo stars. 
Among them are the classical analysis of solar silicon Unes by 
Holweger ("1973^ and the more recent work of Asplund (2000), 
who introduces 3D line formation based on a solar hydrody- 
namical model. 

The introduction of the atomic model used is the topic of 
Section 2. In Section 3 two atmospheric models of the Sun are 
presented. Section 4 shows the results of the NLTE calculations 
obtained with different settings of the remaining free parame- 
ters. In Section 5 spectral line synthesis under both LTE and 
NLTE assumption is applied to determine the photospheric so- 
lar Si abundance and the Si i/Si ii ionisation equilibrium. The 
results are discussed in the final section. 

2. Atomic model of silicon 

As usual, the level system of silicon was extracted from the 
NIST data banlfl. Silicon in cool stars is fully represented by 
the first three ionisation stages, where for the purpose of NLTE 
calculations Sim reduces to the ground state, 3s^ 'S. The con- 
figuration of the neutral atom includes 4 electrons, with that 
of the lowest energy being 3s'^3p^. This gives rise to singlets 
and triplets with a different outer electron replacing 3p. A sin- 



gle quintet term is found near 4 eV, which is due to a 3s 3p^ 
configuration, and which does not couple to any other singlet 
or triplet term with known line strength, except for the 3/7 ^P 
- 3p^ ^S°. This transition is so weak that the quintet term can 
safely be neglected. A similar configuration exists in ionised 
silicon, where the quartet term 3s3p^'^P at 5.3 eV is ignored. 
Our silicon atomic system therefore consists only of the sin- 
glet/triplet system of Si i, the doublet system of Si ii, and the 
singlet ground state of Si iii. It is shown in Fig. [1] as an en- 
ergy level diagram with all the terms included here. To keep 
the number of lines and levels manageable, the level system 
was reduced to terms implying a thermal differential distribu- 
tion of the corresponding level populations. Thus, a total of 132 
terms of Sii, 41 terms of Sin plus the Sim ground state give 
our approach a realistic background. As with other elements, 
completeness is lost at high excitation energies. This holds in 
particular for the Si ii ion, where the last eV is not represented. 
However, since these terms are more than 15 eV from the Si n 
ground state, the population of the levels is extremely low. 

Line transitions are based on the work of the Opacity 
Project d 19941 1, in particular the calculations of Nahar & 
Pradhan ( 119931 ). The total number of lines included is 786 for 
Sii and 182 for Sin. For NLTE transfer calculations simple 
Gauss profiles were used with 9 wavelength points each. 

Fortunately, bound-free transition cross-sections with com- 
plex structure are available from calculations at the TOPbas^l 
(see also Nahar & Pradhan '1993). The cross-sections of the 
lowest three terms, 3/?^P, 3/? 'D and 3/?'S are an order of 
magnitude greater than those of the next terms, which is very 
similar to the configuration of Mg (Zhao et al. 119981 ). It tends 
to decouple the metastable terms efficiently from the excited 
ones. 

Background opacities are calculated with an opacity sam- 
pling code based on the hydrogen lines, the line lists made 
available by Kurucz (|1992| l, and on the important bound-free 
cross-sections of hydrogen and the most abundant metals. The 
background opacities are sampled on a random grid of be- 
tween 5000 and 10000 wavelengths, to which are added the 
wavelengths of the line profiles. The final NLTE line formation 
program thus samples roughly 14000 wavelengths. In the UV, 
where most of the background opacity is expected, the wave- 
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Fig. 3. Temperature stratifications of ODF (dotted) and OS 
(solid line) solar models (left axis). Difference between the two 
stratifications, Ar(T) = r(T)OS - r(T)ODF (right axis). 

length intervals of the random grid are between 1 and 4 A. We 
emphasise that the choice of this frequency grid affects only 
the bound-free radiative interaction rates. Earlier experiments 
with a finer frequency mesh of up to 40000 randomly sampled 
frequencies for Fe i (Gehren et al. 1200 11 1 showed that the influ- 
ence of a very dense frequency sampling is less than marginal. 
This is similar to the sampling situation for OS model atmo- 
spheres, where changes in temperature corrections drop below 
a few K when the number of frequencies is increased to values 
beyond 10000 (Grupp'2004, Fig. 8). Background opacities for 
bound-bound transitions are sampled at the corresponding line 
frequencies. 

For our Si NLTE calculations we take into account inelastic 
collisions with electrons and hydrogen atoms for both excita- 
tion and ionisation. The formulae of van Regemorter (119621 ) 
and Allen ( I1973I I are used to describe the excitation of allowed 
and forbidden transitions by electron collisions, respectively. 
Ionisation cross-sections for electron collisions are calculated 
applying the formula of Seaton ( T962b . Drawin's d 19681 H 9691 ) 
formula as described by Steenbock & Holweger ( 119841 ) is used 
to calculate neutral hydrogen collisions, with a similar formula 
for bound-free hydrogen collisions. To allow for some empir- 
ical correction to the Drawin approximations, a scaling factor 
S\i is applied to the formula in our calculations. 

3. Solar model atmospheres 

Solar atmosphere models are supposed to give the basic back- 
ground for the NLTE line formation program as well as to the 
spectrum synthesis. This investigation assumes a standard ID 
model atmosphere in radiative-convective equilibrium, where 
convection is treated in a mixing-length approach. Since there 
exist two different opacity approximations based either on 
opacity distribution functions (ODF) or on opacity sampling 
(OS), it seems worthwhile extending our research to the differ- 
ences of the silicon level populations between the two atmo- 
spheric models. 



For this purpose both solar models are considered. The 
ODF model is our old MAFAGS standard atmosphere, nearly 
identical to that of Kurucz ( |1992| l. but downscaled from 
Kurucz' value of logepe = 7.67 to a solar Fe abundance of 
logspe = 7.51. It was used in all previous NLTE calculations 
(e.g. Gehren et al. '2001). The OS model is based on a sim- 
ilar model atmosphere code of Grupp ( 120041 ). but with con- 
vection treated according to Canuto Sl Mazzitelli ( 119921) . and 
with opacities sampled from continuous opacities and from 
an extended line list of Kurucz & Bell ( 119951 ) at 86000 ran- 
domly chosen frequency points. The two types of atmospheres 
lead to slightly different temperature stratifications, as can be 
seen in Fig. [3] The two models differ in temperature outside 
logT ^ 0.3 with a moderately increasing amplitude. Note that 
the differences outside logr - -3 do not affect the observed 
Si spectrum, whereas inside log t - 0.3 the gap between the 
temperature stratifications is the result of different versions of 
the mixing-length theory (see Grupp, 2004, for a more detailed 
comparison). The latter marginally affects the formation of the 
Balmer lines, but does not influence Si line formation very 
much. 

4. NLTE line formation results 

The Si departure coefficients for the more important terms of 
Si I and Si ii shown for our final selection of the hydrogen col- 
lision parameter (Sh = 0.1) in Fig. |4] are surprising in that they 
follow LTE out into the middle photosphere. This has not been 
encountered in any of the other neutral ions analysed so far 
It is the result of an unusually large energy gap between the 
ground state, 3p^P, and the two metastable terms 3/5 'D, and 
3/9 'S, on one side, and the first excited levels of neutral sili- 
con on the other. These gaps are ~ 5, 4, and 3 eV, respectively, 
and they shift all lines emerging from those levels into the UV, 
most of them below 2000 A, where the Al i absorption edge 
is suppressed the radiation field by an order of magnitude. As 
a result, all these strong lines have only weak radiative rates 
due to a very low integrated mean radiation field. There is a 
marginal similarity with the Mg i ion (Zhao et al. I1998I I. how- 
ever, with the difference that the more important low-excitation 
transitions in neutral magnesium are well above 2000 A and 
thus considerably stronger This is also the case for Fe i which 
shares ionisation energy and solar abundance with both Mg and 
Si but has a much more complex line spectrum with strong lines 
in the visible. 

Therefore, the interaction of the three most populated Si i 
levels with the Si ii ion is completely based on bound-free pro- 
cesses, where photoionisation and ionisation by electron colli- 
sions compete in strength, with a small net pumping to the Si ii 
ground state. As would be expected, the 3p 'D and 3/7 'S terms 
follow 3/9 tightly. In the absence of stronger radiative rates 
the excited terms are determined by collisions in their respec- 
tive regions of line formation. 

Test calculations applying increased photoionisation or de- 
creased electronic collisions reveal that Sii is relatively sta- 
ble against a variation of atomic data. Although a strongly in- 

^ according to their population and /-value 
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Fig. 4. Departure coefficients bi - Ni/Ni^^^ for the final 
parametrization of the Si atomic model (Sh = 0.1) calculated 
for an ODF (top) and for an OS (bottom) solar model atmo- 
sphere as function of the standard optical depth. Note the sim- 
ilarity of the results for both atmospheric models. The doublet 
terms 4s ^S, 4p belong to Si ii. 



Fig. 5. Departure coeflicients for the final atomic model in the 
ODF solar atmosphere. All parameters are the same as in Fig. 
m except that the photoionisation rates were scaled by a factor 
of 10 (top), and the electronic collision rates divided by a factor 
of 10 (bottom). 



creased photoionisation (see Fig.|5l top) starts to depopulate the 
levels already in deeper layers, this does not affect neutral Si 
line formation very much, because most of the lines observed 
in the visible and near infrared solar spectrum are weak and 
thus formed in deeper layers. The only strong line at 3905 A 
retains a thermalised source function out to the upper photo- 
sphere. The ionisation itself changes only by a very small frac- 
tion. The bottom of Fig.|5]demonstrates a variation of the elec- 
tron collision rates by a factor of 0. 1 . Similar to the increased 
photoionisation, though less pronounced in its influence on the 
departure coefficients, the reduction of the electronic collision 
efficiency also leads only to a marginal de-thermalisation, that 
does not affect neutral silicon line formation very much. Only 
the near-infrared lines respond to collisions (see below). 

While the Sin 3p^P° ground state is a true mirror of the 
Si I ground state, the excited levels of Si ii lead to absorp- 
tion lines coming from much deeper layers in the atmosphere. 
This holds in particular for the highly excited 4s ^S - 4p^P° 
transition with lines at 6347.095 and 6371.360 A, respectively. 
Although absorbed at an excitation energy of more than 8 eV, 
the lines are clearly visible in the spectra of solar-type stars and 



may therefore serve as a test for the ionisation and Si ii excita- 
tion equilibrium in the Sun. 

A variation of the hydrogen collision scaling factor 5^ be- 
tween and 1 allows us to introduce similar changes as those 
resulting from electron collisions. Due to the large number of 
hydrogen atoms, the neutral hydrogen collision rates can be 
significantly stronger than the electron collision rates. Yet, in 
the range of 5^ parameters considered here, they do not lead 
to significant changes in the line profiles in the visible where, 
generally, abundance differences for the two extremes (Sh- 
and 1) stay within 0.03 . . . 0.04 dex. This situation changes sig- 
nificantly, whenever the Sii near-infrared lines are analysed. 
As is evident from Table 1 most of these are absorbed from 
either 4s ■'P", 4p 'P, or 4p^D. Since all the excited terms are 
only loosely coupled to each other and to other high-excitation 
terms, their departure coefficients tend to diverge from each 
other at logr ^ -2, i.e. where the cores of the relatively 
strong near-infrared lines form. This is the reason for a slightly 
increased abundance trend going with the hydrogen collision 
scaling factor For lines beyond 1 fim the corresponding abun- 
dance variations can reach more than 0. 1 dex, and their cores 
suggest the choice of 5h = 0. 1 . 
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All calculations were carried out with a revised version of 
the DETAIL program (Butler & Giddings 1985 ) using acceler- 
ated lambda iteration following the extremely efficient method 
described by Rybicki & Hummer ( [19911 [r992j . 

5. Line synthesis and the solar Si abundance 

The observed solar flux spectrum was taken from the Kitt Peak 
Atlas (Kurucz et al. |1984| l. Spectral line synthesis is performed 
interactively using the SIU program package of Reetz (1991). 
The profiles of all the lines given in Table 1 were calculated 
for various assumptions, with atomic parameters varied as de- 
scribed above. Van der Waals damping constants log Ce for Si i 
lines were computed according to the Anstee & O'Mara ( 19911 
11995 ) interpolation tables. For use in other stars the log^/e© 
values in the last column are most important. They are given in 
Table 1 for our final atomic model statistical equilibrium cal- 
culations with Sh = 0.1, and they are used to determine solar 
Si abundances below. Van der Waals broadening constants for 
the Si II lines are purely empirical allowing an optimal fit of 
the line profiles. Since the profiles are simultaneously affected 
by macroturbulence broadening, similar fits could be obtained 
also with ~ 5 times smaller damping constants. Due to this bias 
the resulting abundances are uncertain by 0.05 dex. 

Fig.|6]shows some fits of the near-infrared lines between 1 .0 
and 1 .3 fj.m. In contrast to the Si lines in the visible, the cores of 
the near-infrared lines deviate from the LTE profiles. It is only 
this variation with the degree of thermalisation that allows us 
to determine a best fit value of 5^ = 0. 1 . As mentioned above, 
this selection does not affect the visible line profiles. 

Nearly all of the Si line profiles observed in the solar spec- 
trum display some kinematic asymmetry with a red-shifted line 
core compared to the fits shown in Fig.|6]. This difference does 
not depend on line blends, for some of which information is 
evidently missing. The effect is accompanied in most lines by a 
corresponding blue shift of the blue line wings, and sometimes 
by observed line wings that are both deeper than any profile 
fit possible under line synthesis with the currently adopted so- 
lar parameters. We note that a large part of such defects are 
the result of our approximate handling of photospheric veloc- 
ity fields. In a ID analysis of spectral line formation, all known 
atmospheric motions are replaced by a micro-/macroturbulence 
approach that can at best account for a mean velocity field. A 
3D analysis would probably remove most of these shortcom- 
ings (see Asplund 2005. for a review). Modelling the complex 
granular velocities with a ID model, however, would at least re- 
quire the macroturbulence velocities H to be adjusted to the re- 
gions of line formation, which are different for weak and strong 
lines. This is impossible with a single macroturbulence convo- 
lution, because higher S values would be needed for the weak 
lines than for the strong line wings. 

When it comes to spectral line synthesis, every effort to de- 
termine the absolute solar photospheric silicon abundance is 
hampered by uncertain or inconsistent /-values. Since the life- 
time experiments of Savage & Lawrence (I1966I I and Marek & 
Richter 09731) there has been an ongoing discussion of the 
measurements regarding the lifetime of the Sii 45 level. 
While all other level lifetimes led to consistent results, this 
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Fig. 7. Solar abundances calculated with three different sets of 
/-values. The mean value is shown as a straight line only for 
the Garz (1973) data, together with the dotted +lcr error lines 



showed controversial values depending on the experimental 
method used (Bashkin et al. 1980; Becker et al.T980ai Becker, 
Zimmermann & Holweger il980bi Smith et al. il987j . The life- 
time measurements of the 4^ ^Pj level were carried out with 
three methods: phase shift, beam foil, and laser fluorescence, 
of which the last, used by Becker et al. (I1980al) . should be the 
most reliable. However, this value, r = 4.4 ± 0.4 ns is at vari- 
ance with the four independent lifetime results obtained by the 
other authors, 6.8 + 0.8, 5.95 + 0.48, 5.9 + 0.7 and 6.3 + 0.6 
ns, respectively. Thus the correspondingly calibrated absolute 
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Fig. 6. Near-infrared Sii solar line profiles. NLTE and LTE profiles are calculated for the same abundance, where the NLTE 
profile with = 0. 1 refers to the best fit. 



/-values differ by roughly 25%, which translates into a system- 
atic abundance difference of 0.10 dex. This is still the current 
systematic uncertainty of the photospheric solar sihcon abun- 
dance. 

To determine the solar Si abundance, it is also necessary 
to judge the available /-values according to their consistency. 
In Fig. |2] both the Kurucz and the NIST /-values produce 
an enormous abundance scatter with standard deviations of 
cr(log Sq) - 0.40 and 0.38, respectively. On the other hand, 
the relative /-values of Garz ( |1973l l are still the only source for 



high-quality data of lines in the visible, with log eo(Garz) - 
7.62 ± 0.05. None of the sets of oscillator strengths displays a 
clear trend in either of the different parameter ranges, although 
it may be significant that the three most deviating abundances, 
based on Kurucz & Peytremann's oscillator strengths, emerge 
from the 4s^P° term, i.e. the term for which the lifetime mea- 
surements were most contradictory. Removing five outlyers of 
the NIST abundances from the sample, the resulting scatter 
about the mean abundance is significantly reduced, now with 
log eo(NIST) =7.52 + 0.13. 
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Fig. 8. Solar profile fits for highly excited Si n lines at 6347 
and 6371 A. Blend calculations do not include terrestrial atmo- 
spheric lines. 

As seen in Table 1, it is possible in the solar spectrum to 
measure the line profiles absorbed from the Si ii 4i^Si/2 level, 
which is excited to more than 8 eV. The results are shown in 
Fig. m the line data and Si abundances are found in Table 1 
and Fig. |7] respectively. Both LTE and NLTE give the same 
abundances within 0.01 dex. Although not perfect, the ioni- 
sation equilibrium, i.e. the abundance difference between Si ii 
and Si i (using Garz' /-values) is ^ -0.08 dex, just at the limit 
of a 2cr interval. It is important to note that the two lines de- 
pend sensitively on the choice of the damping constant, repre- 
sented here by C^. Using Kurucz' van der Waals broadening 
constant which is more than a factor 10 smaller, and replac- 
ing the opacity difference by slightly stronger blends, implies a 
change in the abundances by ^ 0.05 ... 0. 10 dex. Thus the Si ii 
abundances add to the results of the neutral data, but they are 
by no means decisive. 

Altogether, applying the laser fluorescence lifetime mea- 
surements to the Garz /-values would increase the latter by 0. 1 
dex and thus reduce the abundance by the same amount. The 
resulting common mean abundance together with the NIST val- 
ues would be log Sq = 7.52 + 0.09. This would include the Si ii 
lines and lead to a perfect ionisation equilibrium. 

The influence of NLTE and of the choice of the model at- 
mosphere on the solar abundances is best viewed using the line 



Fig. 9. Logarithmic solar abundance differences for silicon 
lines. 

wavelength as an x-axis. Here we present in Fig. the log- 
arithmic abundance differences between NLTE (final atomic 
model, Sn = 0. 1) and LTE, both from profile fits using the ODF 
solar model. Although the mean abundance difference is small, 
the systematic increase towards near-infrared wavelengths is 
remarkable. The lines with low NLTE abundances all have rela- 
tively strong equivalent widths between 100 and 350 mA. Their 
/-values are much stronger than those of the visual lines that 
arise from the same lower levels, mostly 4-s^P°, 45 'P°, and 
4p ^D, because the visual lines lead to even higher excited lev- 
els (see Fig. While most of the visible lines are formed 
well inside log rgooo = -2, the region of line formation for the 
near-infrared lines is shifted outward to log rgooo = -3 ... - 2, 
where the line source function of those lines differs signifi- 
cantly from the local Planck function. For the strongest near- 
infrared lines, however, the core contributes only marginally to 
the abundance determination, which is dominated by the strong 
line profile wings coming from a line formation region inside 
logTgooo = -2. Therefore, the trend in Fig.|9t documents the 
different line saturation (see also Fig.|6]l. 

Fig. |9]3 displays the abundance differences calculated for 
the final NLTE models between the ODF and the OS solar at- 
mosphere. Here, we note a similar trend with wavelength, but 
the corresponding correlation with the strong lines is much less 
pronounced. It is more likely that the different photospheric 



J.R. Shi et al.: Statistical equilibrium of silicon in the solar atmosphere 



9 



7.70 r 



S 7.60 r 




7.40 r • 



7.30 r H 

I t \ t t t \ t t t \ t t t \ t t t \ t 

0.4 0.6 0.8 1.0 1.2 

Wavelength [|lm] 

Fig. 10. Corrected logarithmic solar Si abundances with /- 
values from Garz (119731 filled circles, corrected by +0.1) or 
from NIST (filled triangles, omitting eight outlyers). See text 
for discussion. 

temperature gradients (see Fig. O contribute to the change in 
the abundance differences, because the weaker Si i lines in the 
visible are formed at 50 K higher temperatures and thus require 
a slightly higher abundance. The mean Si abundance does not 
depend on the choice of the model atmosphere. 

6. Discussion 

Both the atomic model and the line formation results do not 
differ substantially from the previous analysis of Wedemeyer 
( 1200 11 1. Our aims of atomic model completeness were more 
directed towards inclusion of high energy levels, whereas 
Wedemeyer apparently favoured level completeness up to a 
limiting energy. Therefore, he included a Si i quintet and some 
Si II quartet terms, all of which are coupled only by weak for- 
bidden transitions to the leading multiplets. In fact, our test cal- 
culations show that these terms are insignificant for the depar- 
ture coefficients. 

Therefore it is not too surprising that our results agree in 
many details. Most important is our agreement with respect to 
the nearly thermal population of all Sii levels inside optical 
depth log T5000 = -2. In fact, Wedemeyer' s departure coeffi- 
cients (see his Fig. 2) are even more thermalised than ours. The 
remaining differences probably come from slightly different in- 
teraction cross-sections. Our solar silicon abundances are also 
very similar, except for an offset that is mostly due to his ad- 
justment of the Garz (119731) oscillator strengths to the Becker 
et al. (ll980a|[T980bl) lifetime measurement. We also agree with 
his result that in the solar photosphere Si abundance corrections 
due to NLTE departures are small. However, it must be noted 
that such abundance corrections become significant in the near 
infrared, where the cores of strong lines are important. 

As is the case for all elements, the agreement between me- 
teoritic abundances and those obtained from photospheric anal- 
ysis is a major topic of research. Because the solar photospheric 
silicon abundance is used as a reference for ff-elements in the 
Galaxy, and as the normalising element for meteoritic abun- 



dances in the solar system, it is even more important for the 
early chemical evolution of the Milky Way. In particular the 
question arises as to whether the meteoritic and photospheric 
Si/Fe ratios are the same. Fig. [10] combines the selected re- 
moval of eight lines with deviating NIST /-values with the 
correction of the Garz /-values for a 25% smaller lifetime of 
the Si I 45 ^^Pj level as measured by Becker et al. ( 1980a). This 
brings together the blue, visual and near-infrared wavelength 
region on a common abundance scale and avoids excessive 
scatter The resulting photospheric abundance, loge0(cor) - 
7.52 + 0.06 is in good agreement with the results of Asplund 
(I2OOOI) and Wedemeyer ( 200TT l. 

During the last years it has become evident that ID abun- 
dance analyses carry systematic errors as compared with 3D 
results. While the 3D results are only a few, and their reliability 
is not yet beyond doubt, it seems that some of the metals anal- 
ysed in 3D or 2D lead to slightly different photospheric abun- 
dances. One of the main obstacles against a full hydrodynamic 
investigation is still the limits of computing speed and memory, 
respectively. Yet preliminary calculations cited by Wedemeyer 
and the 2D simulations of Steffen & Holweger ((2002i) indi- 
cate solar abundance corrections for weak Si i lines of ^ -1-0.02, 
whereas Asplund (120001 1 favors a 3D hydrodynamical correc- 
tion of ^ -0.02 for his set of Si lines, both with reference to the 
solar model of Holweger & Miiller ( 119741 ). The final abundance 
thus may be between 7.50 and 7.54. With a meteoritic ratio of 
logesi - logspe = 0.06 (Lodders'2003', Asplund et al. l2005T l. 
the meteoritic Fe abundance would be as low as 7.44 ... 7.48. 
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Table .1. Atomic data for spectral lines of Si i and Si ii (last two entries). References to the ^/-values are G: Garz ( 119731 
experimental), K: Kurucz & Peytremann ( 119751 ). and N: NIST ( |http://www.physics.nist.gov/| ). Values for log gfe^ refer to our 
final atomic model, with 5h = 0.1. 



Mult 


A 


Line 


^low 


logCfi 




log 




[A] 






[eV] 










2 


4102.93 


ip Sq 


A ^ 3 DO 

— 45 r J 


1.909 


-30.972 


1. i A 

— :>.14 


(j 


4.52 


3 


3905.53 


ip So 


A ^ loo 

— 45 r ] 


1.909 


-30.917 


1 no 
—L.yjy 


(j 


6.41 


4 


12103.56 


A ^ 3r>o 


A „ 3r\ 

-4p D, 


4.707 


-30.759 


-U.Jj 




7.15 




11991.57 


A r. 3 DO 




A r. 3TA 

-4p D, 


4.920 


-30.760 


—yj.Zy 


IN 


7.26 




12031.54 


A ^ 3 DO 


y1 . . 3r\ 

-4p D3 


4.954 


-30.743 


r\ A ^ 
U.41 


JN 


7.93 




12270.72 


A ^ 3 DO 

4^ 


- 4p U2 


4.954 


-30.752 




XT 


7.13 




12395.85 


A r. 3 DO 

4^ 


-4p D, 


4.954 


-30.757 


—L.JJ 


XT 

JN 


5.83 


5 


10603.45 


A ^ 3 DO 

45 


/I .. 3d 

-4p 


4.707 


-30.677 


—U.J 2 


XT 

JN 


7.16 




10749.40 


A ^ 3 DO 

45 


-4p P, 


4.707 


-30.689 


— U. j4 


XT 

JN 


7.30 




10786.88 


A r. 3 DO 

45 r ^ 


A ^ 3o 

-4p Pq 


4.707 


-30.691 


— U.4Z 


XT 

JN 


7.16 




10661.00 


A ^ 3 DO 


y1 .. 3d 

-4p P, 


4.920 


-30.687 


r\ A ^ 
—U.41 


XT 

JN 


7.22 




10827.10 


A c 3 DO 

45 r-, 


A 3o 

-4p P, 


4.954 


-30.677 


U. 1 / 


XT 

JN 


7.71 




10979.34 


A r. 3 DO 

45 


/I ^ 3d 

-4p P, 


4.954 


-30.688 


— U. 


XT 

JN 


6.95 


6 


10371.30 


A ^ 3 DO 

45 


-4p S, 


4.707 


-30.659 


f\ A ^ 


XT 

JN 


6.65 




10288.90 


A r. 3 DO 

45 


/I r. 3 C 

-4p S, 


4.920 


-30.661 


—yj.oy 


XT 

JN 


5.85 




10585.17 


A r. 3 DO 

45 


-4p Sj 


4.954 


-30.659 


— U. io 


XT 

JN 


7.36 


9 


5793.07 


A r. 3 DO 

45 r ^ 


-5p 


4.707 


-30.294 


—z.Uo 




5.57 


10 


5645.61 


A c 3 DO 

45 r ^ 


^ 3o 

-5p Pj 


4.707 


-30.294 


O 1 A 

— z. 14 





5.47 




5690.43 


A c 3 DO 

45 r ^ 


^ 3d 

-5p P, 


4.707 


-30.294 


— 1.0 / 





5.77 




5701.11 


A f 3po 
45 r ^ 


n 30 

-5p Pq 


4.707 


-30.294 


— Z.Uj 





5.55 




5665.55 


A r. 3 DO 

45 


^ . . 3o 

-5p P, 


4.920 


-30.294 


o r\A 
— Z.U4 




5.51 




5708.40 


A r. 3 DO 

45 


^ 3d 

-5p Pj 


4.954 


-30.294 


1 AH 

—1.4/ 




6.12 


11 


5622.23 


A f 3 DO 
45 r ^ 


n 3 C 


4.707 


-30.294 


— 1.D4 


XT 
JN 


4.46 




5684.48 


A c 3 DO 

45 r-, 


^ 3 c 

-5p S, 


4.954 


-30.294 


— L.OJ 




5.93 


11.04 


4818.03 


A c 3 DO 

45 r-, 


3r^ 

-6p D3 


4.954 


-30.294 


— L.J 1 


XT 

JN 


4.93 


16 


5948.54 


A r. loo 

45 r ^ 


-5p D2 


5.082 


-30.629 


— l.Zo 


yj 


6.43 


17 


5772.15 


A r. loo 

45 r ^ 


-5p Sq 


5.082 


-30.287 


— 1. Ij 




5.89 


17.09 


4947.61 


A r. iOO 

45 r ^ 


-6p Sq 


5.082 


-30.287 


101 
—l.ol 


XT 

JN 


5.26 


22 


7415.96 


ip 


"4/ P3 


5.616 


-29.969 


— U. jU 


K 


6.80 


23 


7405.79 


ip D, 


-4/ Pj 


5.614 


-29.869 


— U.oZ 




6.88 


23 


7415.36 


ip 


/I /■ 3t7 

-4/ P2 


5.616 


-29.869 


— 1 . oU 


K 


5.84 


26 


7226.21 


ip D, 


-4/ i->2 


5.614 


-29.800 


1 ^ 1 

— 1. jl 




6.09 


27 


6244.47 


1 r.3 3tao 

ip 


-5/ l->2 


5.616 


-29.868 


— u.oy 




6.22 


28 


6237.32 


1 r.3 3r^o 

ip D, 


^ /■ 3t7 

-5/ Pj 


5.614 


-29.869 


—U.JO 




6.43 


28 


6243.82 


1 ^3 3r^o 

ip 


-5/ P3 


5.616 


-29.868 


-U. / / 


K 


6.26 


29 


6145.02 


1 ^3 3r^o 

ip 


-5/ <-l3 


5.616 


-29.869 


—U.oZ 


K 


6.12 


29 


6155.14 


ip D3 


^ /■ 3/^ 


5.619 


-29.869 


— U.4U 


Jv 


6.73 


30 


6125.02 


1 r.3 3"r\o 

ip D, 


^ /■ 3r^ 


5.614 


-29.869 


n m 


K 


5.98 


30 


6142.49 


1 r.3 3"r\o 

ip D3 


^ /■ 3r^ 


5.619 


-29.869 


— u.yz 




6.03 


31 


10843.87 


4p'Pi 


-4rf'D° 


5.863 


-30.145 


0.15 


N 


7.42 


32 


10627.66 


4p'P, 


-4rf3p° 


5.863 


-30.692 


0.00 


K 


7.11 


36 


7680.27 


4p'P, 


-5d'D° 


5.863 


-29.656 


-0.69 


G 


6.94 


43 


8752.01 


3^/'D° 


-4/'F3" 


5.871 


-30.751 


-0.52 


K 


7.19 


44 


8742.46 


3rf'D° 


-4/3F3 


5.871 


-30.916 


-0.63 


K 


7.10 


50 


7034.91 


3rf'D° 


- 5/^63 


5.871 


-29.027 


-0.88 


G 


6.63 


53 


10689.73 


4p3D, 


-4rf3F° 


5.954 


-29.964 


0.01 


N 


7.42 


53 


10694.27 


4p3D, 


-4d^F° 


5.964 


-29.944 


0.16 


N 


7.56 


53 


10727.43 


4p3D3 




5.984 


-29.907 


0.27 


N 


7.75 


53 


10784.57 


4p^^D, 




5.964 


-29.965 


-0.72 


N 


6.81 


53 


10882.83 


4p3D3 


_4rf3po 


5.984 


-29.945 


-0.73 


N 


6.84 


57 


7918.38 


4p3D, 


-5rf^F° 


5.954 


-29.663 


-0.61 


G 


6.98 


57 


7932.35 


4p3D2 




5.964 


-29.663 


-0.47 


G 


7.15 


57 


7944.00 


4p3D3 


-5dX 


5.984 


-29.663 


-0.31 


G 


7.30 


57 


7970.30 


4p^^D, 


-5d^F° 


5.964 


-29.663 


-1.47 


G 


6.19 



J.R. Shi et al.: Statistical equilibrium of silicon in the solar atmosphere, Online Material p 3 



Table 1 (conlinucd) 
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